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We study the effect of absorptive corrections due to parton recombination on the 
parton distributions of the proton. A more precise version of the GLRMQ equations, 
which account for non-linear corrections to DGLAP evolution, is derived. An analysis of 



Abstract 



HERA Fi data shows that the small- x gluon distribution is enhanced at low scales when 



the absorptive effects are included, such that a negative gluon distribution at 1 GeV is no 
longer required. 

x. 

g3 1 Introduction 

At very small values of x it is expected that the number density of partons within the proton 
becomes so large that they begin to recombine with each other. This phenomenon of parton 
recombination is also referred to as absorptive corrections, non-linear effects, screening, shad- 
owing, or unitarity corrections, all leading to saturation. The first perturbative QCD (pQCD) 
calculations describing the fusion of two Pomeron ladders into one were made by Gribov- 
Levin-Ryskin (GLR) [1] and by Mueller-Qiu (MQ) [2]. The GLRMQ equations add an extra 
non-linear term, quadratic in the gluon density, to the usual DGLAP equations for the gluon 
and sea-quark evolution. The evolution of the gluon distribution is then given by 

dxg(x,Q 2 ) _ a s , 9 a|(Q 2 ) f l dx' { 2 ]2 

a'=q,g ^ 

where R ~ 1 fm is of the order of the proton radius. The GLRMQ equations account for all 'fan' 
diagrams, that is, all possible 2 — > 1 ladder recombinations, in the double leading logarithmic 
approximation (DLLA) which resums all powers of the parameter asln(l/x) \u(Q 2 /Qq). 



There has been much recent theoretical activity in deriving (and studying) more precise 
non-linear evolution equations, such as the Balitsky-Kovchegov (BK) and Jalilian-Marian- 
Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equations (see [3] for a review). Note 
that the BK and JIMWLK equations are both based on BFKL evolution. However, for the 
most relevant studies in the HERA and LHC domain (x > 10~ 4 ), the predominant theoretical 
framework is collinear factorisation with DGLAP-evolved parton distribution functions (PDFs). 
At very small values of x it might be expected that the DGLAP approximation would break 
down, since large ash^l/x) (BFKL) terms would appear in the perturbation series in addition 
to the «5 ln(<5 2 /Qo) terms resummed by DGLAP evolution. However, it turns out that the 
resummed NLL BFKL calculations of the gluon splitting function P gg [4] and the gluon trans- 
verse momentum distribution [5] are rather close to the DGLAP calculations. Moreover, the 
convolution P gg ® g(x,Q 2 ) coincides with the NNLO DGLAP result and is close to the NLO 
DGLAP result for x > 10~ 4 [6]. Hence, in the analysis of current data, it is reasonable to ignore 
BFKL effects. 

If recombination effects are significant, it is therefore important that they be incorporated 
into the global DGLAP parton analyses which determine the PDFs from deep-inelastic scatter- 
ing (DIS) and related hard-scattering data. Such a programme, based on GLRMQ evolution 
(which accounts for gluon-induced screening only), was implemented some years ago [7], before 
the advent of HERA. The input gluon and sea-quark distributions were assumed to have a 
small-x behaviour of the form xg,xS ~ x~ 0,5 at an input scale of Ql = 4 GeV 2 . The inclusion 
of shadowing effects, both in the form of the input PDFs and in the GLRMQ evolution, was 
found to significantly decrease the size of the small- x gluon distribution in comparison with 
the result with no absorptive corrections. A crucial observation is that, at that time (1990), 
F 2 data were only available for x B > 0.07, and so these results were largely dependent on the 
theoretical assumptions made for the starting distributions. However, with HERA, we now 
have F 2 data down to x B ~ 10~ 4 or less, and so the PDFs at small x can be determined directly 
from the HERA data. 

In fact, the advent of HERA data has led to a puzzling behaviour of the small-x gluon 
and sea-quark PDFs at low scales Q 2 . If we write xg ~ x~ Xg and xS ~ x~ Xs , then the 
expectation of Regge theory is that X g = \ s = X soit for low scales Q < Q ~ 1 GeV, where 
A so ft — 0.08 [8] is the power of s obtained from fitting soft hadron data. At higher Q > 1 
GeV, QCD evolution should take over, increasing the powers X g and As. However, the current 
MRST2004 NLO [9] and CTEQ6.1M [10] PDF sets exhibit a very different behaviour at low 
scales from that theoretically expected; see Fig. [TJ In fact, the MRST group has found that a 
negative input gluon distribution at Qo = 1 GeV is required in all their NLO DGLAP fits since 
MRST2001 [11]. The CTEQ group, who take a slightly higher input scale of Q = 1.3 GeV, 
also find a negative gluon distribution when evolving backwards to 1 GeV. 

Since data at small x B now exist, the introduction of the absorptive corrections is expected 
to increase the size of the input gluon distribution at small x to maintain a satisfactory fit to 
the data. To understand this, note that the negative non-linear term in the GLRMQ equation 

slows down the evolution. Therefore, it is necessary to start with a larger small- x gluon 
distribution at low scales Q ~ Q to achieve the same PDFs at larger scales required to describe 
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Figure 1: The behaviour of the gluon and sea-quark distributions at Q 2 = 2 GeV 2 found in the 
MRST2004 NLO [9] and CTEQ6.1M [10] global analyses. The valence-like behaviour of the 
gluon is evident. 

the data. If the non-linear term is neglected, the input small- a; gluon distribution is forced to 
be artificially small in order to mimic the neglected screening corrections. 

We have anticipated that the introduction of absorptive corrections will enhance 1 the small- 
x gluon at low scales, and hence could possibly avoid what appears to be anomalous behaviour 
at small x. Thus, here, we perform such a study using an abridged version of the MRST2001 
NLO analysis [11], improving on our previous analysis [13]. First, we derive a more precise 
form of the GLRMQ equations. 

2 Non-linear evolution from diffractive DIS 

The inclusive proton structure function, F 2 (x B) Q 2 ), as measured by experiment, can be approx- 
imately written as a sum of the single Pomeron exchange (DGLAP) contribution and absorptive 
corrections due to a 2 — > 1 Pomeron merging; see Fig. |2fa)- That is, 

F 2 (x B ,Q 2 ) = F 2 DGLAP (x s ,Q 2 ) + AFt(x B ,Q 2 ). (2) 

In computing AF 2 abs we need to sum over all possible cuts. The Abramovsky-Gribov-Kancheli 
(AGK) cutting rules [14] were originally formulated in Reggeon field theory but have been 
shown to also hold in pQCD [15]. Application of the AGK rules gives the result that relative 
contributions of +1, —4, and +2 are obtained according to whether neither Pomeron, one 
Pomeron, or both Pomerons are cut; see Fig. |2fb). Therefore, the sum over cuts is equal to 
minus the diffractive cut and so the absorptive corrections can be computed from a calculation 

1 Eskola et al. [12] have found that taking input gluon and sea-quark distributions at Q 2 = 1.4 GeV 2 , then 
evolving upwards with the GLRMQ equations based on LO DGLAP evolution, improves the agreement with 
F 2 data at small x B and low Q 2 compared to the standard CTEQ sets, and leads to an enhanced small- a; gluon 
distribution for Q 2 < 10 GeV 2 . Note, however, that there is a large NLO correction to the splitting function 
P qg which changes completely the relationship between the quark and gluon distributions, and so weakens the 
conclusion of Ref. [12]. 
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Figure 2: (a) Absorptive corrections to F 2 due to the 2 — > 1 Pomeron contribution, (b) 
Application of the AGK cutting rules. For simplicity, the upper parton ladder, shown in the 
right-hand diagram of (a), is hidden inside the upper blob in each diagram of (b). 

of the t-integrated diffractive structure function F 2 D ^ 3 \x]p, (3, Q 2 ), where (3 = x B jxjp and xjp is 
the fraction of the proton's momentum transferred through the rapidity gap. 

The pQCD description of F 2 D ^ is described in [16,17]. Working in the fixed flavour number 
scheme (FFNS), it can be written as 



,D(3) 



? D(3) 
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+ p m) + F 
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apart from the secondary Reggeon contribution. The separation between the soft Pomeron and 
QCD Pomeron is provided by a scale /x ~ 1 GeV. For simplicity, we take /x to be the same 
as the scale Qo at which the input PDFs are taken in the analysis of F 2 data, so //q = Qq = 1 
GeV, the value used in the MRST2001 NLO analysis [11]. The contribution to the absorptive 
corrections arising from the soft Pomeron contribution of (J3J) is already included in the input 
PDFs, therefore 
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where / p .di S s. is the fraction of diffractive events in which the proton dissociates. In practice, 
we take / p .di S s. — 0.5 and take an upper limit of 0.1 instead of 1 for xjp in (jlj). 2 

First consider the contribution to (J3J) from the F^^ rt term. 3 It corresponds to a 2 — > 1 



2 The value of / p .diss. = 0.5 is justified by a ZEUS comparison [18] of proton-tagged diffractive DIS data with 
data which allowed proton dissociation up to masses of 6 GeV, where /p.diss. = 0-46 ± 0.11 was obtained. 
3 The other two contributions to (@J) are described after (| 1 3D . 
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Pomeron merging with a cut between the two Pomeron ladders and can be written as 



pert. ) 



(5) 



where C^a are the same coefficient functions as in inclusive DIS. The diffractive PDFs, aP = zq° 
or zg B , where z = x/xjp, satisfy an inhomogeneous evolution equation [17]: 



a pert.( x iP, z , Q 2 ) = ! A* 2 ) fliP ( z ' ^ 



ua pert. 

dlnQ 2 



2tt 



,/ g) a 



/D 
pert. 



Paip(z) flp(Xjp] Q 2 



(6) 
(7) 



a'=g,g 



Here, fjp(xjp\ Q 2 ) is the perturbative Pomeron flux factor, 
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(8) 



The diffractive slope parameter Bp, comes from the t-integration, while the factor R g accounts 
for the skewedness of the proton gluon distribution [19]. There are similar contributions from 
(light) sea quarks, where g in (JHJ) is replaced by S = 2(u + d + s), together with an interference 
term. A sum over all three contributions is implied in (JB|) and in the second term of (|7J). The 
Pomeron PDFs in (JBJ), a IP {z 1 Q 2 ] fi 2 ), are evolved using NLO DGLAP from a starting scale 
fi 2 up to Q 2 , taking the input distributions to be LO Pomeron-to-parton splitting functions, 
<F{z^ 2 ^ 2 )=P oIP (z) [17]. 
From 



a(x, Q z 



.DGLAP, 



x,Q 2 ) + Aa abs (a;,Q 



where a(x, Q 2 ) = xg(x,Q 2 ) or xS(x, Q 2 ), and 

1 Z" 1 

Aa abs (x,g 2 ) = -— — / da;^ a% ert X x iPi x / x JPiQ d 

/p. diss. Jx 



(9) 



(10) 



Differentiating @ with respect to Q 2 gives the evolution equations for the (inclusive) gluon 
and sea-quark PDFs: 



da(x, Q 2 
dlnQ 2 



^ V P aa ,®a ! 



o =3,9 



l-/p. 



p. diss. Jx 



dxp Pajpi^x/xp) fjp(xjp; Q 2 ). 



(11) 



Thus (jUj) is a more precise version of the GLRMQ equations ((T}, which goes beyond the 
DLLA and accounts for sea-quark recombination as well as gluon recombination. Consider 
the recombination of gluons into gluons, for example, in the DLLA where x <C Xjp, then 
Pgip = 9/16 [17]. Taking R g = 1 and / p .di S s. = 0, then (fTTj) becomes 



dxg(x, Q 2 
dlnQ 2 



2tt 



P 9 a ' 
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9 a 2 s (Q 2 ) pdxjp 
16 B D Q 2 



Xjp 



[xjpg(xjp, Q' 



2M 2 



(12) 
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Comparing to (JTJ) this is simply the GLRMQ equation with R 2 = 8B^- For numerical results we 
take Bd = 6 (4) GeV -2 for light (charm) quarks, which would correspond to R = \Z8Bd = 1.4 
(1.1) fm. 

The procedure for incorporating absorptive corrections into a (NLO) global parton analysis 
(in the FFNS) is as follows: 

1. Parameterise the x dependence of the input PDFs at a scale Qo ~ 1 GeV. 

2. Evolve the PDFs xg(x,Q 2 ) and xS(x,Q 2 ) using the non-linear evolution equation (|11)1. 
(The non-singlet distributions are evolved using the usual linear DGLAP equations.) 

3. Compute 

F 2 (x B , Q 2 ) = ® a ~ 737 f dx *> feet + <S.4l , (13) 

a=qg 1 Jp.diss. A fl 

and compare to data. Here, the two terms inside the square brackets are beyond collinear 
factorisation, that is, they cannot be written as a convolution of coefficient functions 
with the PDFs. The first term inside the square brackets corresponds to the process 
7*ZP — > cc. The second term corresponds to the process 7*/P — > qq, for light quarks 
with a longitudinally polarised photon. These contributions are calculated as described 
in Ref. [17]. 

As usual, these three steps should be repeated with the parameters of the input PDFs adjusted 
until an optimal fit is obtained. This procedure is our recommended way of accounting for 
absorptive corrections in a global parton analysis. However, in practice, available NLO DGLAP 
evolution codes, such as the qcdnum [20] program, are often regarded as a 'black box', and it 
is not trivial to modify the usual linear DGLAP evolution to the non-linear evolution of (fTTj) . 
Therefore, we adopt an alternative iterative procedure which avoids the explicit implementation 
of non-linear evolution, but which is equivalent to the above procedure. 



3 Effect of absorptive corrections on inclusive PDFs 

We model our analysis of HERA F 2 data [21] on the MRST2001 NLO analysis [11], which 
was the first in which a negative gluon distribution was required at the input scale of Qo — 1 
GeV. (The more recent MRST sets have not changed substantially at small x.) We apply cuts 
x B < 0.01, Q 2 > 2 GeV 2 , and W 2 > 12.5 GeV 2 , leaving 280 data points. The input gluon and 
sea-quark distributions are taken to be 

xg(x,Q 2 ) = A g x- X3 (l-x) 3 - 70 (l + e g Vx~ + -y g x) - x' 5 ~ (1 - x) 10 , (14) 
xS(x,Q 2 ) = A s x' Xs (l-x) 7 - 10 (l + e s V^ + lsx), (15) 

where the powers of the (1 — x) factors are taken from [11], together with the valence-quark 
distributions, uy and dy, and A = d — u. The A g parameter is fixed by the momentum sum 
rule, while the other nine parameters are allowed to go free. Since we do not fit to DIS data 
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with x B > 0.01, we constrain the input gluon and sea-quark distributions, and their derivatives 
with respect to x, to agree with the MRST2001 NLO parton set [11] at x — 0.2. This is done 
by including the value of these MRST PDFs at x — 0.2, and their derivatives, as data points 
in the fit, with an error of 10% on both the value of the MRST PDFs and their derivatives. 
Therefore, the PDFs we obtain are not precisely constrained at large x, but this paper is 
primarily concerned with the small- x behaviour of the PDFs. 
The procedure we adopt is as follows: 

(i) Start by performing a standard NLO DGLAP fit to F 2 data with no absorptive corrections. 

(ii) Tabulate AF 2 abs , given by (jlj), and Aa abs , given by (JTUJ), using PDFs g(xjp, fi 2 ) and 
S(x]p,fj, 2 ) obtained from the previous fit. 

(iii) Perform a standard NLO DGLAP fit to 'corrected' data, F 2 DGLAP = F 2 -AF 2 abs , to obtain 
PDFs a DGLAP . Then correct these PDFs to obtain a = a DGLAP + Aa abs . These latter 
PDFs a then satisfy the non-linear evolution equations (jllj) . 

(iv) Go to (ii). 

Each successive iteration of steps (ii) and (iii) introduces another level of 2 — > 1 Pomeron 
mergings, so that eventually all the 'fan' diagrams are included, achieving the same effect as 
the procedure described at the end of Section |21 

Note that the correction to the PDFs, a = a DGLAP + Aa abs , in each step (iii), was omitted in 
our previous analysis [13]. Consequently, the effect of the absorptive corrections on the PDFs 
at large scales was overestimated. Also in [13], the known LO Paip(z) were multiplied by free 
parameters ('K-factors'), determined from separate fits to diffractive DIS data, in an attempt 
to account for higher-order pQCD corrections to the LO Pomeron-to-parton splitting functions. 
However, since these K-factors took unreasonable values, with some going to zero, here we have 
chosen to fix them to 1. Therefore, the updated analysis, presented here, does not require a 
simultaneous fit to the diffractive DIS data. 

In Fig. Efa) we show the gluon distribution at scales Q 2 = 1, 4, 10, and 40 GeV 2 obtained 
from fits before and after absorptive corrections have been included. Both fits are almost equally 
good with x 2 /d.o.f. values of 0.86 and 0.87 for the fits without and with absorptive corrections 
respectively. At low Q 2 the absorptive corrections give an increased gluon distribution at 
small x, apart from at x < 10 -4 where there are only a few data points and where additional 
absorptive effects (Pomeron loops) may become important. The non-linear term of (fTTj) slows 
down the evolution, so that by 40 GeV 2 the two gluon distributions are roughly equal; see 
Fig. Ha). 

We repeated the fits without the negative term in the input gluon distribution, that is, 
without the second term in (fT4^) . When absorptive corrections were included, almost the same 
quality of fit was obtained (x 2 /d.o.f. = 0.90), while without absorptive corrections the fit was 
slightly worse (x 2 /d.o.f. = 0.95). We conclude that absorptive corrections lessen the need for 
a negative gluon distribution at Q 2 = 1 GeV 2 . The gluon distributions obtained from six 
successive iterations of steps (ii) and (iii) above are shown in Fig. E^b). The convergence is 
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Figure 3: (a) The gluon distribution obtained from fits to F 2 data, before and after absorptive 
corrections have been included, (b) The effect of successive iterations on the gluon distribu- 
tion obtained from fits to F 2 , taking a positive definite input gluon at 1 GeV. Each iteration 
introduces another level of 2 — > 1 Pomeron mergings. 



fairly rapid, with only the first three iterations having a significant effect, that is, the 'fan' 
diagrams which include 8 — > 4 —>■ 2 — > 1 Pomeron mergings. 

Although we have seen that the inclusion of absorptive corrections has reduced the need for 
a negative gluon, it has not solved the problem of the valence-like gluon. That is, the gluon 
distribution at low scales still decreases with decreasing x, whereas from Regge theory it is 
expected to behave as xg ~ x ~ Xsott with A so f t ~ 0.08. We have studied several possibilities 
of obtaining a satisfactory fit with this behaviour [13]. The only modification which appears 
consistent with the data (and with the desired X g = As equality) is the inclusion of power-like 
corrections, specifically, a global shift in all scales by about 1 GeV 2 . (Note that a similar 
shift in the scale is required in the dipole saturation model [22].) However, we do not have a 
solid theoretical justification for this shift. Therefore, a more detailed, and more theoretically- 
motivated, investigation of the effect of power corrections in DIS is called for. 
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